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The neutron elastic magnetic form factor was extracted from quasielastic electron scattering on 
deuterium over the range — 1.0 GeV^ —4.8 GeV^ with the CLAS detector at Jefferson Lab. High 
precision was achieved with a ratio technique and a simultaneous in-situ calibration of the neutron 
detection efficiency. Neutrons were detected with electromagnetic calorimeters and time-of-flight 
scintillators at two beam energies. The dipole parameterization gives a good description of the 
data. 

PACS numbers: 14.20.Dh, 13.40.Gp 



The elastic electromagnetic form factors of the pro- 
ton and neutron are fundamental quantities related to 
their spatial charge and current distributions. The dom- 
inant features of the larger form factors G^j, G^, and 
G^Ij were established in the 1960's: the dipole form 
Gd = (1 + QVA)-2 where A = 0.71 GeV^ gave a good 
description of these form factors {G^j/fip « G^^/fin ~ 
Gd) within experimental uncertainties, corre- 
sponding (at least for ^ 1 GeV'^ or large radii) to 
an exponential falloff in the spatial densities of charge 
and magnetization. Recent Jefferson Lab results on the 
proton form factors show a dramatic departure from the 
dipole form even at moderate fi\ while the neutron 
magnetic form factor falls below the dipole at high 



Q2 {Glj/nnGo = 0.62 ± 0.15 at = 10 GeV^ 0). 
Describing all these modern results with nucleon models 
and lattice calculations has been a challenge [3,0)S[1)0]- 
Also, the elastic form factors are the zeroth moment of 
the generalized parton distributions (GPDs) and thus, 
constrain GPD models Last, we note that some 

models predict significant deviations from the dipole for 
Q2 <5 GeV2 1,1:^. 

To distinguish among different models, high precision 
and large coverage are important. At larger momen- 
tum transfer G^j is known much more poorly than the 
proton form factors 8] . In this Letter we report on a new 
measurement of G^j in the range = 1.0 — 4.8 GeV^ 
at Jefferson Lab. The precision and coverage of these 
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results eclipse the world's data in this range. System- 
atic uncertainties were held to 2.5% or less. 

In the absence of a free neutron target, we measure 
the ratio R of the cross sections for the ^H(e, e'ri)p and 
^H(e, e'p)n reactions in quasielastic (QE) scattering on 
deuterium. A nucleon with most of the momentum 
from the scattered electron is detected in coincidence 
with the final state electron. The ratio R is defined as 
R = ^[2H(e,e'n)Q£;]/^pH(e,e'p)Q£] & K [H [11 
and 
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and separated by 4.7 cm. The downstream cell was filled 
with liquid hydrogen for calibrations and efficiency mea- 
surements. The upstream cell was filled with liquid deu- 
terium for the ratio measurement. The target was made 
of aluminum with 20-micron aluminum windows. The 
CLAS vertex resolution of 2 mm enabled us to separate 
events from the different targets [l^ . 

We now describe the analysis. Nucleons from 
quasielastic events tend to be ejected close to the di- 
rection of the 3-momentum transfer q while inelastically 
scattered nucleons are not We required the an- 

gle 0pq between the nucleon 3-momentum and q to be 
small (61™""^ = 2.5° - 4.5° across the range) and 
(1) integrated over all azimuthal angles about q. Another 



where E is the beam energy, aMott is the cross sec- 
tion for scattering off a scalar (spinless), point parti- 
cle of unit charge, r = Q^/AM'^, M is the nucleon 
mass, and 9 is the electron scattering angle. The fac- 
tor a{E, Q^, d^q""^ ,W^g_^) corrects for nuclear effects and 
depends on E and cuts on 0^"^, the maximum angle 
between the nucleon direction and the three-momentum 
transfer g, and W^^^^, the square of the maximum value 
of the mass recoiling against the electron assuming a sta- 
tionary target. We used the one-photon exchange ap- 
proximation in the numerator of Eq. [T] to express the 
cross section in terms of the neutron form factors. The 
right-hand side of Eq. [T] contains the desired G^,j along 
with the better-known proton cross section and the neu- 
tron electric form factor (G^), which is believed to be 
small over the range here. For QE kinematics (within 
a cone 6*™°^ around q) can be extracted from Eq. [1] 
as a function of by relying on knowledge of the proton 
cross section {i.e., the Arrington parameterization [l3|). 
G^, calculations of a{E, Q"^, 0™""^, W^^^), and measure- 
ments of R. The ratio method is less vulnerable to nu- 
clear structure {e.g., choice of deuteron wave function, 
etc.) and experimental effects {e.g., radiative correc- 
tions, etc.). The challenge here is to accurately measure 
the nucleon detection efficiencies. 

The two reactions were measured in the CLAS detec- 
tor [1^ at the same time and from the same target to 



reduce systematic uncertainties. Two electron-beam en- 
ergies were used, 2.6 GeV and 4.2 GeV. CLAS consists 
of six independent magnetic spectrometers each instru- 
mented with drift chambers I5|, time-of- flight (TOF) 
scintillators covering polar angles 8° < 9 < 143° [16|, a 
gas-filled threshold Cherenkov counter (CC) and a 
lead-scintillator sandwich-type electromagnetic calorime- 
ter (EC) covering 8° < < 45° 0. CLAS was triggered 
on electrons by requiring a coincidence between CC and 
EC signals in one sector. Neutrons were measured sepa- 
rately in the TOF and EC. Protons were measured using 
the drift chambers and TOF systems. A novel dual-cell 
target was used consisting of two coUinear cells each 5- 
cm long - one filled with -'^H and the other with - 



cut, < 



1.2 GeV eliminated most in- 
cut. Our simu- 



elastic events that survived the 0^°^ 
lations of the quasielastic [§| and inelastic production 
[lit show the fraction of inelastic events surviving these 
cuts is less than 0.5% of the total. To measure R accu- 
rately, the sohd angles of CLAS for the ^H(e, e'njQE and 
2H(e, e'p)QE reactions have to be identical. The nucleon 
solid angles were matched by first determining event-by- 
event the nucleon momentum from the electron kinemat- 
ics assuming quasielastic scattering. The expected pro- 
ton and neutron trajectories in CLAS were checked to 
see if both trajectories would lie within the CLAS accep- 
tance. Only events where both nucleons were expected 
to strike CLAS were analyzed. 

Once the event sample was selected, corrections for 
the detector efficiencies and other effects were applied. 
Neutrons were measured in two CLAS scintillator-based 
detectors: the EC and the TOF. The neutron detec- 
tion efficiency (NDE) measurement was performed using 
tagged neutrons from the ^Il{e,e'Tr'^)n reaction, where 
the mass of the unobserved neutron was inferred from 
the measured electron and pion kinematics and matched 
with possible hits in the neutron detector. The value 
of the detection efficiency can vary with time-dependent 
and rate-dependent quantities like photomultiplier tube 
gain so it was measured simultaneously with the primary 
deuterium measurement. The measured neutron detec- 
tion efficiency for each sector for the TOF and for nine 
subsections in each EC sector were fitted with polyno- 
mials at low neutron momenta and a constant at high 
momenta. The EC efficiency typically reached a maxi- 
mum value of RiO.6 while the maximum TOF efficiency 
was ftiO.OS [^, HO]. The calibration target was also used 
to measure the proton detection efficiency using elastic 
scattering p{e, e'p). The kinematics of the scattered elec- 
tron were used to predict the location of the elastically 
scattered proton in CLAS and the event was searched for 
a proton at that location. 

The calculation of the nuclear correction factor, 
a(S,Q2,0™"^, W^„^), in Eq.fflis described in Ref. [2l|. 
The cross section was calculated using the Plane Wave 
Impulse Approximation (PWIA) for > 1.0 GcV^, the 
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FIG. 1: (color online). Results for 7? as a function of for 
two beam energies. Each set is a weighted sum of the TOF 
and EC neutron measurements. Only statistical uncertain- 
ties are shown. Numerical results are reported in the CLAS 
Physics Data Base [20l ]. 



AV18 deuteron wave function [22], and Glauber theory 
for final-state interactions (FSI). The correction is the ra- 
tio of the full calculation to the PWIA without FSI. The 
correction was averaged over the same 9pq range used in 
the analysis and was less than 0.1% across the full 
range. 

In our analysis we assumed QE kinematics and ignored 
the Fermi motion that can knock the ejected nucleon out 
of the acceptance. To correct for this effect we simulated 
QE scattering from a fixed target nucleon and tested to 
see if it struck the active area of CLAS (an 'expected' 
event). We then simulated the nucleon's internal mo- 
tion (with the Hulthen distribution) and elastic scatter- 
ing from this moving particle. With the target momen- 
tum known (in simulation) we re-calculated the trajec- 
tory to see if it still struck CLAS and satisfied the 6'™°^ 
cut (an 'actual' event). The ratio of actual to expected 
events is the correction for that nucleon. The ratio of 
these corrections for the neutron and the proton multi- 
plies R. The correction to is in the range « 0.9 — 1.3. 

We present our results for R in Fig. [1] for the two beam 
energies and for > \ GeV^ where we have overlap- 
ping TOF and EC data. The corrections described above 
have been included and only statistical uncertainties are 
shown. For each beam energy we averaged the two neu- 
tron measurements (EC and TOF) weighted by the sta- 
tistical uncertainties. Measurements of R at the same 
but different beam energies are not expected to be 
the same because the kinematics are not the same (re- 
call Eq. [1]). The data cover the range with excellent 
statistical accuracy and with a large overlap between the 
two data sets. 

A detailed study of each correction's contribution 
to the systematic uncertainty has been performed [J]. 
Listed in Table H] are the largest contributions along with 
the maximum (typical) value across the full range. 



The largest contributions come from the parameteriza- 
tions of the neutron detection efficiencies for the TOF 
and EC systems. To estimate the uncertainty associated 
with the NDE measurement, the order of the polynomial 
and position of the edge of the constant region used to fit 
the data were varied to determine the effect on G\j as a 



function of Q 

The extraction of G\[ 



Uncertainties were in the range 0.5-3.2%. 

depends on the other elastic 
form factors (see Eq.[T|) and their uncertainties contribute 
to the uncertainty in G^j. The proton cross section un- 
certainty was estimated using the difference between pa- 
rameterizations by Bosted and Arrington 13, [i^]. The 
average difference was < 1% with a maximum of 1.5%. 
For G^, the difference between the Galster parameteri- 
zation and a fit by Lomon was used 2J, |25| with a max- 
imum uncertainty of 0.7%. The upper limit of the 9pq 
cut was varied by ±10%, changing G^ by a maximum 
of about 1.0% and by 0.3% on average [9|. The uncer- 
tainty of the Fermi motion correction was calculated us- 
ing two dramatically different momentum distributions 
of the deuteron: a flat one and the Hulthen distribution. 
This correction to G^j changes by < 1% between the two 
Fermi motion distributions. The quadrature sum of the 
remaining, maximum systematic uncertainties was less 
than 0.5% The final systematic uncertainty for the 
EC measurement was < 2.4% and for the TOF measure- 
ment it was < 3.6%. 

The CLAS extraction of G\j{Q^) consists of overlap- 
ping measurements. The TOF scintillators cover the full 
angular range of CLAS, while the EC system covers a 
subset of these angles, so G^j can be obtained from two 
independent measurements of the e — n production. The 
experiment was performed with two beam energies with 
overlapping coverage so the detection of nucleons of 
a given occurs in two different regions of CLAS. Four 
measurements of G^ have been obtained from CLAS 
that could have four semi-independent sets of system- 
atic uncertainties. Shown in Fig. [2] are the results for 
G^^ from the different measurements divided by /i„GD 
for normalization and to reduce the dominant depen- 
dence. Only statistical uncertainties are shown. Here 
the different measurements should agree because G\j 
depends only on Q^. The two measurements for each 
beam energy are consistent within the statistical uncer- 
tainties, suggesting the systematic uncertainties are well- 
controlled and small. The results in Fig. [2] were then 



Quantity 




Quantity 




EC NDE 


< 1.5% (1%) 


TOF NDE 


< 3.2% (2%) 


proton (7 


< 1.5% (0.8%) 




< 0.7% (0.5%) 


Fermi loss 


< 0.9% (0.5%) 


6pq cut 


< 1.0% (0.3%) 


Remainder 


< 0.5% (0.2%) 







TABLE I: Upper limits (typical values) of systematic errors. 
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combined in a weighted average as a function of . The 
final systematic uncertainty varied from 1.7-2.5% across 
the full data range. The larger uncertainty on the pa- 
rameterization of the TOF NDE (see Table 1) did not 
push the total, weighted uncertainty above our goal of 
3%. There are more calorimeter data due to its higher 
efficiency and the maximum EC uncertainty was 1.5% 

The final, combined results for G^j are shown in Fig. [3] 
with a sample of existing data [3,[ll[2i,t27, 28, 29]. The 
uncertainties are statistical only. Systematic uncertain- 
ties are represented by the band below the data. A few 
features are noteworthy. First, the quality and cover- 
age of the data is a dramatic improvement of the world's 
data set. Second, our results are consistent with previous 
data, but with much smaller uncertainties. Third, the 
dipole form is a good representation here, which differs 
from parameterizations and some calculations at higher 
where previous results for G^lj/ {^nGo) decrease with 
increasing [1, 0, 0] ■ We note there appears to be an 
offset between the low-Q^ end of our data and some ear- 
lier results [ll|, [2^ that is about twice the uncertainty of 
the offset. Last, any apparent fluctuations in our results 
{e.g. at 1.29 GeV^) are not significant enough to draw 
any firm conclusions here. 

The curves shown in Fig. [3] are from Diehl et al. 
Guidal et al. [|, and Miller [i] and are all constrained by 
the world's previous data. In Diehl et al. the GPDs are 
parameterized and fitted to the experimental data (green 
band). The curve reproduces some of the low-Q^ data, 
but lies above our results. Guidal et al. use a Regge 
parameterization of the GPDs to characterize the elas- 
tic nucleon form factors at low momentum transfer and 
extend it to higher (dashed line). The curve repro- 
duces the existing, higher data (which fall well below 
the dipole in the range = 6 — 10 GeV^), but is not 
consistent with our results. In Miller's calculation the 
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FIG. 2; (color online). Results for G^i/d-inGo) as a function 
of for four different measurements (two beam energies). 
Only statistical uncertainties are shown. 
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FIG. 3: (color online). Results for Clj / {finG d) from the 
CLAS measurement are compared with a selection of previous 
data [13, [H, m, 113, m. 111] and theoretical calculations [I, i, 
[^. Numerical results are reported in the CLAS Physics Data 
Base ii. 



nucleon is treated using light-front dynamics as a rela- 
tivistic system of three bound quarks and a surrounding 
pion cloud (solid curve). The model gives a good descrip- 
tion of much of the previous data even at high and is 
consistent with our results. 

The neutron magnetic form factor was measured in the 
range = 1.0 - 4.8 GeV^ with the CLAS detector at 
Jefferson Lab using the ratio of e — rt to e — p scattering. 
Two incident beam energies were used and systematic un- 
certainties were < 2.5%. Neutrons were measured with 
two independent systems: time-of-fiight scintillators and 
electromagnetic calorimeters. Detector efficiencies were 
measured simultaneously with the production data using 
a dual-cell target containing and ^H. The data pro- 
vide a significant improvement in precision and coverage 
in this range and are surprisingly consistent with the 
long-established dipole form. The calculation by Miller 
is in good agreement with our results. 
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